Crossover from Luttinger- to Fermi-liquid behavior in strongly anisotropic systems 

in large dimensions 
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We consider the low-energy region of an array of Luttinger liquids coupled by a weak interchain 
hopping. The leading logarithmic divergences can be re-summed to all orders within a self-consistent 
perturbative expansion in the hopping, in the large-dimension limit. The anomalous exponent scales 
to zero below the one-particle crossover temperature. As a consequence, coherent quasiparticles with 
finite weight appear along the whole Fermi surface. Extending the expansion self-consistently to all 
orders turns out to be crucial in order to restore the correct Fermi-liquid behavior. 



PACS numbers : 71.10.Pm, 71.10.Hf, 71.27.+a 

One-dimensional metals are characterized by the fail- 
ure of the quasiparticle concept, signaled by the vanishing 
of the single-particle spectral weight at the Fermi surface 
(FS). This produces a number of additional anomalous 
single- and two-particle properties , as the absence of a 
step in the momentum distribution, spin-charge separa- 
tion, the vanishing of the local density of states p(ui) 
as a power law like p(to) cx o> Q , as well as nontrivial 
power-law behavior of several correlation functions. All 
these anomalous properties (except for spin-charge sepa- 
ration) are controlled by the single exponent a . Metals 
showing these properties are called Luttinger liquids (LL) 
PJ. On the other hand, most isotropic two- and higher- 
dimensional metals show Fermi-liquid (FL) properties, 
their low-energy regime can be described by well-defined 
quasiparticles, and the exponents of the correlation func- 
tions coincide with those inferred from dimensional anal- 
ysis. 

We thus encounter two distinct universality classes 
when going from one to higher dimensions, and it's cer- 
tainly a fascinating issue to study the crossover region be- 
tween them. Castellani, Di Castro and Metzner B have 
approached this problem by carrying out an analytic con- 
tinuation from one dimension D = 1 to D = 1 + s. They 
have shown that the FL regime is always recovered for 
D > 1, except for broken-symmetry states. 

An alternative way is to approach the crossover for 
the anisotropic case by starting from a D — 1 dimen- 
sional array of uncoupled LL and switching on a small 
single-particle hopping t±/\/~D between them ||. This 
also represents, in principle, the experimental situation, 
for example of the organic quasi-one-dimensional conduc- 
tors. The crucial question is: does the system go over to 
a FL for arbitrarily small t±_ and sufficiently low tem- 
peratures? One expects that, at least for temperatures 
T greater than a certain crossover temperature T\r> the 
system still shows LL properties. In fact, it has been 
shown (see, e. g. Ref. [0), that the region of validity 
of the LL extends down to T\r> t e // = t^^ 01 ' 
Therefore, the interesting, nontrivial region is the one for 
which the characteristic energy scale of the system, E\\ || 



is smaller than t e ff . Only in this region, can a possible 
FL behavior occur. 

Since t± is small, the natural idea is to carry out a 
perturbative calculation in t± . Many results have been 
obtained by restricting to the lowest order or to the first 
few orders in t± [^],f||9| . An expansion in powers of t± has 
been introduced in Ref. IB] and systematized diagram- 
matically by the author |10| . This was done by extending 
a method introduced by Metzner to expand around the 
atomic limit of the Hubbard model Jll| . The lowest order 
in t± for the (inverse-) self-energy has been studied by 
Wen H and by Boies et al. Bj. Within their approxima- 
tion (below referred to as LO), the quasiparticle structure 
is recovered for arbitrarily small t± on most of the FS. 
However, the anomalous exponent a still dominates the 
low-energy (k/t e // <C 1) [0 behavior for small c±_ since 
the self-energy behaves as |k| 1_Q , i. e., its asymptotic be- 
havior is still LL-like. As a physical consequence, close to 
the regions for which c± vanishes, the spectral function is 
incoherent and the quasiparticles are thus poorly defined 
(see Fig. 2). In particular, the quasiparticle weight (or 
residue) Z vanishes as Z oc (t ± c ± ) a /'- 1 ~ a \ However, the 
LO approximation is not sufficient to make definite con- 
clusions not even for small t± . The reason is that the 
perturbation t± is relevant, i. e., each order in t± car- 
ries a strongly divergent power of the energy £|| Q_1 . Any 
perturbative expansion restricted to lowest order is thus 
uncontrolled at low energies £u <C t e ff- This is the reason 
why theoretical results are still contradictory about the 
nature of the ground state in this energy region. Since, 
as discussed above, this is precisely the relevant region 
for a possible FL behavior, it is worthwhile investigating 
it in a controlled way. 

In this Letter, we analyze for the first time unambigu- 
ously the region £\\ <C t e / / by extending the calculation 
of the self- energy to infinite order in t± . Specifically, 
we sum self-consistently the infinite series of loop dia- 
grams of Fig la. The sum of these diagrams with fully 
dressed hopping (cf. Fig. lb) gives the exact result in 
the (anisotropic) large-D limit. Drastic changes are ob- 
tained with respect to the LO result, which corresponds 
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to taking just the first term on the r.h.s. in Fig. la 
[BfoJ- Our main result (Fig. 2) is that the anomalous 
exponent a scales to zero and the asymptotic behavior of 
the correlation functions is dominated by mean-field-like 
exponents and thus it becomes FL like. In particular, 
our result restores coherent quasiparticles with a finite 
weight Z around the regions where cj_ = 0. This strik- 
ing outcome points out the importance of extending the 
expansion self-consistently to all order in the cumulants. 
Stopping the series at any finite order would still get an 
anomalous behavior and a vanishing Z. 



4 V 



4 f 



I 

+ 



= 00 



(a) 

(b) 
(c) 



FIG.l. Diagrams contributing to the inverse-self-energy T(k) 
(gray circle) in the D — oo limit within an expansion in the 
dressed hopping T± (a), (b) Dressed hopping (dashed line) 
and its diagrammatic expression in terms of the bare hopping 
t± (full line), (c) Dyson equation for the full Green's function 
ff(k) (cf. Ref. 0). 

We consider an hypercubic array of LL with the total 
dimension D, characterized by an exponent a, and with 
a hopping tx/VD coupling them ||. Here, we take the 
spinlcss case since it allows for crucial simplifications in 
the calculation. Since we are interested in the effects and 
in the fate of the anomalous exponent a, we believe that 
spin-charge separation should not play a role. 

Since it is not possible to sum all diagrams in the ex- 
pansion, we want to select a workable subset of diagrams 
according to some physical limit in order to avoid an ar- 
bitrary choice. The D — > oo limit allows us to restrict 
ourselves to the diagrams indicated in Fig. la. Our 
D — > oo procedure is different from the standard dy- 
namical mean-field theory pl| in several aspects. First, 
and most important, in our case the system is strongly 
anisotropic, since the hopping in one (in the |j) direction 
is not rescaled by the \fD factor and is much larger than 
in the other D — 1 (_L) directions. Our system thus repre- 
sents a 1-D chain embedded in an effective self-consistent 
medium. Since this captures the one-dimensionality and 
the anomalous exponent of the initial system, it is par- 
ticularly suited to study the problem of the crossover 
from one to higher dimensions. Notice that the system 
has more degrees of freedom than the standard D — > oo 
single-impurity problem. As a consequence, the self- 
energy is local with respect to the _L coordinates but 
has a nontrivial dependence on the || ones B. The sec- 
ond difference is merely procedural: instead of solving 
the dynamical mean-field equations, we sum all the di- 



agrams contributing to the self-energy in D — > oo. The 
internal lines of these diagrams must be evaluated self- 
consistently (see Fig. 1). Third, since we are interested 
in the asymptotic behavior, we only keep the leading sin- 
gularities in each diagrams, as we will discuss below. 

As anticipated, each power of t± carries a divergent 
power of the energy oc @70- A crucial point in 

our approach is to consider an expansion in the dressed 
hopping 71 instead of the bare one t± (see below and 
Fig. 1). Since the scaling behavior of T± turns out to 
be Tl oc this cancels exactly the power-law diver- 

gences met in the t± expansion. However, logarithmic di- 
vergences, arising from the integration of the cumulants 
still occur, as we shall see below. 

Further ingredients of the perturbation theory are the 
bare (i. e., tj_ = 0) _L-local m+l-particle cumulants of the 
LL II £ c °(x • • •x m |xo' • • -x TO '). For these, we take the 
form valid for T — and large distances Jll| . It will prove 
convenient to work in real (i. e., x) space, in contrast to 
usual perturbation theory. As discussed above, in the 
D — > oo limit, the inverse-self-energy F(xo) is _L-local, 
and is obtained as the sum of the loop diagrams in Fig. 
la as 

r(x ) = s°(xoio) + ]T (-ir / [n d2 y« 

m=l Jl i- m k=l 

xd 2 x fc TL(-x fe ,0) 0°(yo + x , ■ • • ,y m +x m |y , • • • ,y m ) , 



where "1 J, m" indicates the region |xi| > |x2 1 > ■ • ■ > 
|x m | to which the integral can be restricted by exploit- 
ing the symmetry under exchange of the coordinates 
l,---,m. In Eq. (0), Tj_(x, x± — 0) is the dressed 
hopping written in real space. As anticipated, we are 
interested in the dominant low-energy behavior (here, 
l x o|^e// S> 1) of correlation functions and thus we can 
restrict to the leading logarithmic divergences in the loop 
integrals (|l|) juj . These can be shown to come from fur- 
ther restricting to the region |x p | < |xo| for each p, and 
|Xp| < \y q + £ i x « - Yq> - e 2* q > I for each p>q,q',q^ q\ 
and €i = 0, 1, which will be referred to as "0 JJ. m". For 
convenience, we introduce the "restricted renormalized 
cumulants" (RRC) in this region defined as 

) (2) 
= G° (yo + xo, • • • , Ym + x m |y , • • • , y m ) 

- / d 2 x m+ id 2 y m+ iTL(-x m+ i,0) 

J0i}.m+1 

xGc(yo + x , • ■ • ,y m+ i +x m+ i|y , • ■ ■ ,y m +i) ■ 

Comparing (Q) and @), one can verify that T(xo) = 
&(xo|0). 

Within the leading-logarithmic approach discussed 
above, one can show that the RRC can be written in 
terms of the bare cumulants as 

GciYo + ) (3) 

= F m (l , l m )Gc(yo + x , • • • , y m + x m |y , • • • , y m ) 



2 



where the renormalization factors T m are functions of 
the relative coordinates x, through ij = alog(|xi|t e //) 
p3| . We can thus first carry out the integration over 
the "center-of-mass" coordinate y m +i in (0) by simply 
considering the effect on the the bare cumulant. The 
leading logarithmic contributions to this integral can be 
evaluated in a lengthy but straightforward way |l3| ] . 

We now need to evaluate the dressed hopping 71 to 
insert in (H). This can be obtained in k space in terms 
of r(k) pi] by solving the Dyson-like equation of Fig. lb 

1 

7I(k,cx)=txcx/(l-r(k)i x c x ). (4) 

Within our self-consistent calculation |L5[ one must use 
the fully dressed inverse-self-energy ([I]) in (|J). We 
are interested in the leading contribution to 71 (x, 0) at 
large |x|. This comes from two types of singularities in 
momentum space: (i) from the shifted poles of (Q) at 
c± 7^ 0, |k| 7^ and (ii) from the c± = 0, |k| = power- 
law singularity of T(k). Taking, to start with, the LO 
form for T, namely T(k) « £°(k) - |k| Q -\ the shifted- 
pole contribution (i) turns out to behave like |x| 2q ~ 4 , 
while (ii) yields the strongest contribution oc |x|"~ 3 . This 
is due to the fact that the shifted-pole singularity (i) is 
smeared out after integration over the _L momenta. 

We now evaluate 71 in terms of the inverse self-energy 
r, so that it can be inserted in (||). From (||) with 
to = 0, r(x ) = £°(x |0).Fo(7o)- By taking the FT of 
this expression into momentum space, inserting in (jj), 
and transforming 71 again into real space, we obtain, 
within our leading- logarithm scheme 

Tx(xi,o) = ^| Xl r-V ar s** [Mu)+rd(M] , (5) 

where F (h) = 1/Fa(!i), and T (k) = jj-T {k). Thanks 
to the |x| a_3 term in (|), the FT of 71 now scales as 
|k| 1_Q , as already anticipated. Therefore, the diagrams 
written in terms of the dresses propagator 71 all scale in 
the same way for low energies and thus only logarithmic 
singularities are left |lq| . 

Using Eqs. (|2|), and (||), and integrating over 

Ym+1 and x m+1 , yields a recursive self-consistent equa- 
tion [ |l3| for the renormalization functions J- m , which we 
have solved together with the self-consistency condition 

(I) = 1/J-q(1), by expanding in powers of the variable I 
up to 40 th order and analyzed the series by means of sev- 
eral types of differential Pade approximants. The result 
of this analysis yields for the large-positive-/ behavior of 
the renormalization function Toil) ~ e cl , where the ex- 
ponent c turns out to be equal to 1 within a great degree 
of accuracy (about 0.01%). 

Introducing this result in the expression for the inverse- 
self-energy (Eq. (||) with to = 0) yields T(x) — > 
C?°(x|0)(|x| t e ff) a — > l/|x|, i. e. the anomalous ex- 
ponent a exactly cancels out!. The same things hap- 
pens in momentum space. Here, one obtains T(k) = 



^ c °(k)^ (alog |* e ///k|) -» {iLa-rk^)- 1 (since £ c °(k) oc 
|k| a /(iw — rfcii)). This can be interpreted as the LL ex- 
ponent a renormalizing to zero for energies smaller than 

W/ ■ 

The Green's function of the coupled array of LL 
is given by the Dyson-like equation |l(| £(k) = 
(r(k) _1 — t±c±) 1 (cf. Fig. lc). The main impor- 
tant consequence of our result is that the low-energy 
behavior of correlation functions, as, say, the Green's 
function, is now FL like. Specifically, by an homoge- 
neous rescaling of all energies k — * sk, and c± — > scj_, 
t/(k) oc for s 1, i. e., the anomalous dimen- 
sion vanishes. This result reflects onto the spectral prop- 
erties in the region close to ci = 0. The quasi- 




FIG. 2. Quasiparticle weight Z (a) as a function of the off- 
chain kinetic energy c± (in units of t± ) for the LL exponent 
a = 1/3. Spectral function A(u>) (b) for ku = and c± = 0.1, 
and imaginary part of the self-energy ImT~ (c) as a function 
of w. Our D — oo results (solid lines) re compared with the 
LO result (dashed). In (a) we also plot for comparison the 
result obtained including just the first two diagrams on the 
r.h.s. of Fig. la (dotted). 

particle weight Z(cj_) on the FS [&|| F (cx)] is given by 

Z{c±y 1 = ■^l/g(k^ F (c ± ),iuj) iuJ _o, and is plotted in 
Fig. 2a for a = 1/3 and compared with two results ob- 
tained within the LO approximation (dashed line) and by 
taking in addition the first loop diagrams of Fig. la (dot- 
ted), respectively. For small c± , the LO result would give 
a Z vanishing as Z(c±) oc t e ff a c± a ^ 1 ~ a \ thus yielding 
poorly defined quasiparticles around c±_ = 0. Inclusion 
of the first loop gives a vanishing Z too. On the contrary, 
our result yields a finite Z for c± — > as can also be seen 
from the figure. This also reflects in the spectral function, 
plotted in Fig. 2b. Our result shows a coherent quasi- 
particle peak in contrast to the rather broad and shallow 
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maximum produced by the LO approximation. To under- 
stand this behavior, we have plotted the scattering rate 
(ImT^ 1 ) for fen = in Fig. 2c. As one can see, while 
for the LO approximation ImT" 1 vanishes quite slowly 
(as w 1_Q ), in our result ImT -1 vanishes faster than lin- 
early, as required for well-defined quasiparticles. How- 
ever, one should mention that our calculation, restricted 
to the leading divergences, only yields reliable results for 
ImT^ 1 close to the FS points with c± = 0. In other FS 
points, we cannot reliably state whether ImT^ 1 vanishes 
fast enough or not. Arguments similar to the one of ordi- 
nary perturbation theory jlTj cannot be extended to the 
present case, due to the momentum dependence of the 
vertices in the t± expansion. 

In summary, our result shows that a quasi-one- 
dimensional metal with small perpendicular hopping ij_ 
crosses over from a LL behavior with anomalous expo- 
nent a for energies larger than t e ff and flows to a FL 
fixed point with mean-field like exponents for smaller en- 
ergies in large dimensions. This is due to the LL expo- 
nent a renormalizing to zero for £u <C teff- As a con- 
sequence, well defined quasiparticles are recovered along 
the whole FS. The result contrasts with the LO approxi- 
mation for which the anomalous exponent is still present 
in the asymptotic behavior of the self-energy below t e ff 
and the spectrum is incoherent for small c± . 

In principle, we cannot say whether our result is valid 
also for D = 2 or D = 3. However, it turns out that 
the non-_L-local dressed hopping 7j_(x, x± ^ 0) vanishes 
faster than the _L-local one 7j_(x, 0) for large |x|. Non 
_L-local contributions are thus irrelevant and the present 
result could be thus extended to finite dimensions. This 
is valid for values of a smaller than a critical value ct2 P , 
for which a transition to a two-particle regime takes place 
lipf . In addition, it may possibly exist a critical dimen- 
sion below which non-_L-local diagrams may enhance the 
contribution from the shifted |k|,cj_ =/= poles, and give 
rise to nesting or superconducting instabilities at selected 
regions of the FS. One additional important conclusion 
is that it is crucial to include an infinite number of di- 
agrams in the t± expansion, in order to obtain reliable 
results. 
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